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Nonequilibrium steady state and a strict bound on ballistic transport 
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An explicit matrix product ansatz is presented, in the first two orders in the (weak) coupling 
parameter, for the nonequilibrium steady state of the homogeneous, nearest neighbor Heisenberg 
XXZ spin 1/2 chain driven by Lindblad operators which act only at the edges of the chain. The first 
order of the density operator becomes in the thermodynamic limit an exact pseudolocal conservation 
law and yields - via the Mazur inequality - a rigorous lower bound on the high-temperature spin 
Drude weight. Such a Mazur bound is found a nonvanishing fractal function of the anisotropy 
parameter A for |A| < 1. 
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Introduction. Exactly solvable models which exhibit 
certain generic physical properties are of paramount im- 
portance in theoretical physics, in particular in con- 
densed matter and statistical physics where one of the key 
open issues is the transport in low dimensional strongly 
interacting quantum systems. An example par excellence 
of such models is an anisotropic Heisenberg XXZ spin 1 /2 
chain with a constant nearest neighbor spin interaction 
which, in spite of it being Bethe ansatz solvable [1] , still 
offers many puzzles. For example, at high temperature 
and vanishing external magnetic field, it is not clear even 
if the model exhibits ballistic or diffusive spin transport 
[2] • The question is of long-lasting experimental interest 
[3] . Recently, theoretical study of interacting many-body 
systems has got a new impetus by invoking the methods 
of open quantum systems and master equations ^ in the 
study of quantum transport far from equilibrium [51 [3] . 

We consider the markovian master equation in the 
Lindblad form 

^ = -\[H,p{t)]+Y,'2L,p{t)Ll - {LlLu,p{t)} (1) 
fc 

for an open XXZ spin 1 /2 chain with the Hamiltonian 

H - E(2<-i+i + '^^J^+i + (2) 

where — ± i^^J)' '^^j^^ ~ 1, . . . , n are Pauli oper- 
ators on a tensor product space (C^)**", with symmetric 
Lindblad driving acting on the edges of the chain only 

Li^2 = V(l±M)e/2CT±, i3,4 = V(lT/ik/2^T,f. (3) 

In this Letter we construct an exact nonequilibrium 
steady state (NESS) solution of this model in the regime 
of weak coupling (small e) with the method which seems 
to be un-related to the Bethe ansatz and expresses NESS 
in the form of a homogeneous matrix product opera- 
tor. Such an ansatz has been employed earlier for ex- 
act solutions of classical many-body stochastic processes 



[8j, however with an important distinction that here the 
constructing auxiliary matrices satisfy cubic instead of 
quadratic algebraic relations. Our solution gives birth 
to a spin current related conservation law of an infinite 
chain, which is almost local (pseudolocal) in the metallic 
regime |A| < 1, and hence is used in a Mazur inequality 
[2] to bound the spin Drude weight and prove ballistic 
transport for certain values of anisotropy A, while for 
other we use efficient numerical computation. We observe 
that the graph of the Mazur bound versus A exhibits a 
fractal structure. 

Boundary driven XXZ chain. NESS p^o — hmt_>oo p{t) 
is a fixed point of the flow ([T]) 

-i(adff)poo+ePpoo =0, (4) 
where {a,d H)p := [H,p] and 

P i (1 + ti)V+ + 1(1- m)^- , with (5) 
'D±p := 2afpaf - {af , p} + 2(t^ pa^ - {cr^a^,p}. 

Note the distinct roles of two bath parameters: coupling 
strength e controls the strength of coupling to the spin 
baths at the edges of the chain, while driving strength p 
controls the nonequilibrium forcing due to unequal av- 
erage spin polarizations (or chemical potentials, in the 
spinless fermion formulation [10]) of the two baths. For 
example, in the common derivation of the master eq. ([T]) 
[3] weak coupling (small e) is a standard assumption, 
whereas the linear response physics would be mimicked 
by taking small p at nonsmall e (see e.g. [TTJ[T2]). 

Here we address far-from-cquilibrium physics within 
the regime of weak coupling, so we formally expand NESS 
in terms of the coupling parameter e 

oo 

P^^Y.(ierp(Pl (6) 

Plugging the ansatz ^ into the fixed point condition ^ 
results in an operator valued recurrence relation 

(adi7)p(f) = -Pp(P-i), (7) 



2 



for the sequence {p^^-*}, with the initial condition p*^"^ = 
2~"1, which is an infinite temperature equilibrium state. 

Theorem. Solutions of ([t]) in the first two orders read: 



Z is a non-Hcrmitian matrix product operator 



(si,...,s„)G{ + ,-,0} 



R 



n 



(8) 
(9) 

(10) 



where cr^ = 1 and Ao,A± are triple of near-diagonal 
matrix operators acting on an auxiliary Hilbert space H 
spanned by the orthonormal basis {|L), |R), |1), |2), . . .}: 



An = |L)(L| 



|L)(1| 



|R)(R|+;^cos(rA)|r)(r|, 



oo 



c^sin ( 2 

r=l 

oo 



A |r)(r + l|, (11) 



1 A 



where A = arccos A e MUiM and [x\ is the largest integer 
not larger than x. Constant c e C — {0} is arbitrary, but 
we adopt a choice c = 1 for |A| < 1 (A G M) and c = i for 
|A| > 1 (A e iR) making the matrices ( |Tl| ) always real. 

Proof. First order. From ^ it follows that 2?±1 — 
±2((Tf — crfj, and Vl = 2fi{a\ — a^), hence the first 
order ^ satisfies Q for p = 1 if 



[H, Z] 



(12) 



It is important to observe that the ansatz ( 10 ) does not 
contain any (t| operator, while [H, Z] can only contain 

terms with a single cr|. Let us write H = X]?=i ^j-j+i 



with hj^j+i = 2cr+cr7+i + 2cr7CT++i + Acr|cr|+i. Eq. (12) 



implies that all the terms of [ff, Z] where cr| appear m 
the bulk 1 < j < n should vanish. Looking locally at 
each triple of sites around j this means 



Si S9 S-^ 



i]) 



E tr(a;i,aX^,[Vi,,+/.. 

Sl,S2,S3e{±:0} 

xAsiAs^Asg = 0, for aU ri,r3 e {±,0} 
resulting in eight independent 3-point algebraic relations 



[Ao,A+A_]=0, {Ao,A^} = 2AA+AoA+, 
[Ao,A_A+]=0, {Ao,A2_} = 2AA_AoA_, 
2A[A2,A+] = [A_,A^], 2A[A2,A_] = [A+,Alj, 
2A{A2,A+}-4AoA+Ao 
2A{A2,A_}-4AoA_Ao 



(13) 



{A_, A^} - 2A+A_A 
{A+,Al} - 2A_A+A_. 



Sufficient additional conditions that the boundary terms 
of [H,Z] containing ct^jCT^ result in the exact RHS of 



(12 1, when Z is of the form (10), can be expressed in 



terms of auxiliary boundary vectors (L|, |R) 

(L|A_ = (L|A+A_A+ = (L|A+A^ = 0, 
A+|R) = A_A+A_|R) = A^A_|R) = 0, (14) 
(L|Ao = (L|, Ao|R) = |R). (L|A+A_|R) = 1. 



Verifying (13) and (14), which imply ( [l2| ), for the repre- 
sentation (11) results in trivial trigonometric identities. 

Second order. To prove that ([9| satisfies Q for p = 2 
it is sufficient to show [H, [Z, Z^]] = (X>+ + V^){Z - Z^) 
(a) and [H, [Z-Z'^f] = -(P+-P_)(Z-Zt) (b). These 
relations are implied by [Z, <7\ — cr^] = (2?+ -I- 'D-)Z (a') 
and {Z,al-al} = \{V+-'b_)Z (b'), using the identity 



[H, AB] = [H, A]B + A[H, S], together with which 
we have just proven, Hermitian (anti)symmetrization, 
and the property {V±Zy = V±Z'^ . As 'b± and 
(anti)commutator with u\ — cr^ only act on sites 1 and n, 
and since Z on these sites can o nly contain terms with 



CTj^' and (7,7'", respectively (see (14)), it is sufficient to 
check (a') and (b') for all Z G {1, ct^, cr^, ctj^o-^}, which 
reduce to trivial algebra. QED. 

The theorem constitutes the rigorous part of this Let- 
ter. Its essential ingredient, namely the conservation law 



property (12), shall be applied later to the long stand- 
ing problem of the spin Drude weight. However, for es- 
tablishing the existence of a unique NESS poo several 
mathematical issues still have to be addressed: (i) It is 
s not clear a-priori whether recurrence ([T]) should have 
a unique solution. Suppose uniqueness can be proven 
up to some order p^^"^). Then, the solution of a linear 
equation ([?]) is certainly non-unique up to the addition of 
an arbitrary linear combination of conserved quantities 



Qk, [H,Qk] = 0, namely p 



Y^k^kQk- How- 



ever, as has been checked explicitly by means of computer 
algebra for small n < 6 and any A, a unique set of coeffi- 
cients Uk exists such that Vp'^'P^ is in the image of ad , 
and so the solution of ^ for the next order exists. 
For general n this statement remains a conjecture, (ii) 
For ([9]) to give the complete second order of NESS, one 
thus needs to check in addition that Vp^^^ e Imadi?. 
This has been verified explicitly again only for n < 6, 
and is conjectured to be true for any n. (iii) Conver- 
gence properties of perturbation series ^ are unknown, 
although it should be useful even only as an asymptotic 
series. 

Computation of observables. We note that for any fi- 



nite n, due to their near-diagonality the matrices ( 11 ) can 
be truncated to the d = 2 -I- \ n/2\ dimensional auxiliary 
space Hd spanned by vectors {|L), |R), |1), . . . , |[n/2j)} 
only, reproducing exactly the same matrix product 
operator (10). Let us now describe how to compute 
physical observables in NESS, {A) = trpoo^, to second 
order in e. For example, the spin current from site j to 
j + 1, Jj — i{(J^a~_^_l — cr~a^_^_i) results directly from 



(14), (J, 



i(a+a^.+i 

More interesting are expectations of 
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FIG. 1. Scaled magnetization Mj = {af)/e'^ for strongly 
driven fJ. = 1 XXZ chain of length n = 21 at A = 0.5 (a), 
1 (b), 1.1 (c). Red dots mark analytical results (see text). 



local magnetization and spin-spin correlations which are, 
since Z does not contain any cr^, nonvanishing only in 
the second order in e, hence we define Mj :— (ap/e^ and 
Cj^k ■— (o'fo'fc)/^^- The multiple sums appearing in A/j = 



as 



2-"-itr(T| {//(I + - fi)Z^Z} such 

(a|ZZt) = E^^^,^jL|A,,..A,jRf tr^inLi'^r 
can be evaluated using the transfer matrices (TMs), 
T = (1 - ski) A, ® A„ V = i J2s ® A„ where 
± EE ±1, as e.g. Mj = ^(L|(g)(L|f J-iVf "-J|R)(X)|R), and 
similarly for Cj^k and higher correlations. However, eval- 
uation of such TM products can be drastically simplified 



by observing that, since the matrices (11) have only a 
single element in each row, the subspace of diagonal 
vectors K, spanned by (g) {v), v E {L, R, 1, 2, . . .} is 
invariant under the action of TMs T/C C /C, V/C C K,. 
Using identification (g) |;/) — >■ and defining reduced 
TMs T = T|k:,V = V|k;, reading explicitly 

T =|L)(L| + |R)(R| + 1(|L)(1| + |1)(R|) 



Y,{cos^rX)\r){r\ 



V = 



|L>(1| |1>(R| 



4iJ 

oo 



sm 



1 A 



-1 



A |r)(r+l| 



-l)(r| 



(15) 





r+1 
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(r + l| 



sm 



1 A 



-l)(r| 



(16) 



we obtain efficient TM product expressions 

Mj = ^(L|T^'~iVT"-J», (17) 

Cj- fc = ^2^l|t^-ivt'=-^'"1vt"-'=|r), j < k, 

etc, where, for any finite n, the operators ( 15|16 1 can 
again be truncated to c?-dimensional matrices over Jid- 
Note that with our aforementioned choice of constant c 
the matrix elements of T are always positive. 

Examples. For the isotropic case A = 1, the TMs 
T,V have an effective rank 3, i.e. they close on H3, 



so magnetization profile is - interestingly - found ex- 
plicitly to be linear, with constant long-range correlator 
M, = |(n -t- 1 - 2j), C7,-fc ^ ^t^.J ^ k. Similarly, 
for any A of the form A — cos^nl/m), l,m G Z, i.e. 
for rational A/tt, we find that "Hm+i is invariant under 
T,V, since the off-diagonal chains of ( 15|16 ) are cut, 
i.e. matrix elements at position either (m -t- 1, m -f 2) or 
{m + 2,m + 1) vanish, so T, V can be replaced by to -I- 1 
dimensional matrices of Tj^^^^ , Vj^^^^ , independently 
of the system size n. This allows us to find explicit re- 
sults for some small m, say for A = 1/2 = cos7r/3, Mj — 
f(4(5/8)J-i-4(5/8)"-J-(-8)i-"((-5)J-i-(-5)"-^)). 
For |A| > 1, however, the TMs have always infinite rank. 



with exponentially growing matrix elements, but (17) al- 
low for stable and efficient numerical computations, say 
of Mj in only 0{n^) steps. In Fig. n] we plot magnetiza- 
tion profiles, ranging from fiatlike (for |A| < 1) indica- 
tive of ballistic transport to kink-shaped (for |A| > 1) 
with superexponential (in n) growth of Mj near the baths 
j = l,n indicating the superexponential shrinking of the 
perturbative border for e, and consequently also of the 



M, < 1. Note 



spin current (Jj) = \ne, since |(ct| 
that such kink-shaped profiles Mj have been observed 
numerically for strong driving /i = 1 in conjunction with 
the negative differential conductance p^. One has to 
stress that our result, being perturbative in e, cannot be 
used to address the question of conductivity and spin dif- 
fusion, namely the conductivity {Jj) /Vj{<Jj) oc does 
not exist in the weak coupling limit, i.e. the limits e — 
and /i — )■ do not commute in the thermodynamic limit 
(TL) n — > 00. 

Mazur hound on Drude weight. Nevertheless, our re- 
sult offers very interesting physical application. Namely, 



as it is clear from eq. ( 12 ), the operator Z commutes with 
the XXZ Hamiltonian, apart from the boundary terms. 
For |A| < 1, the matrix product operator (10) makes 
sense even for n = 00, where Z becomes a translation- 
ally invariant pseudolocal conservation law in the spirit 
of Ref. [13] . Since almost all matrix elements (11) are 
smaller than 1 for |A| < 1, it is easy to prove that the 
coefficients of expansion of Z in Pauli spin operator clus- 
ters decay exponentially with the length of the cluster. 

Ref. [Hj used the Mazur inequality [5] to rigorously 
estimate the lower bound on the Drude weight D at high 
temperature [15] in terms of some conservation laws Qk 



lim 

t-^oo 2nt 



dt'{J{t')J) > 



{Qk,Qk) 



(18) 



where J := ^,"=1 Jj^ J{t) ■= e'"*Je-'"\ /3 is the in- 
verse temperature, and {A,B) := 2^"tryl.^^i3 is Hilbcrt- 
Schmidt inner product in the operator space, which de- 
fines the norm ||^|| = a/ {A, A), and conserved quanti- 
ties Qk are chosen to be mutually orthogonal (Qk,Qi) = 
Sk.iWQkW^- However, for XXZ chain all the known lo- 
cal, or 'normalizable', conservation laws [\& are orthog- 
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onal to the spin current (J, Qk 
Nevertheless 



(181 vanishes. 



0, so the RHS of 
many numerical computa- 
tions (e.g. [17]) suggested that for |A| < 1, D should be 
positive, although doubts have been raised recently [5], 
making the issue quite controversial. Here we propose to 
use the first-order-term in NESS (js]), Qz ■= i{Z - Z^) 
as Hermitian operator, which is conserved in TL, since 
clearly hm„_>.oo Qz]||/||Qz|| — 0, and is relevant to 
the current, since (J, Qz) — {n—l)/2. Putting Qz 
and performing TL on both sides p[H] of (18), we obtain 
!^Dz where Dz ■= lim„^oo(»^ - 



D 



lim„ 



Dr, 



> 



lf/[2n{Qz,Qz)] = ilim„^^n/(Z,Z). Computing 
(Z, Z) = 2"" tr ZZt by iterating TM T we get efficiently 
computable expression for the Mazur bound 

1 Th 

Dz = - lim — — . (19) 

4«^oo (L|T"|R) ^ ' 

For lAI > 1, we derive Dz = from the fact that the 



matrix elements in the tridiagonal band of T (151 are 



all positive and exponentially growing, which allows for 
numerically observed spin diffusion at infinite tempera- 
ture [TTJ [12]. Consider now |A| < 1. For a dense set 
of values of the anisotropy A = cos(7r//m), TM T can 
be, as discussed earlier, replaced by a finite m + \ di- 
mensional matrix T' = T|-h„,i, and Dz (19) can be 



evaluated exactly by means of Jordan decomposition of 



T' 



For example, for A 
/l 1/2 Q\ 
10 
1/2 1/4 3/8 
yo 3/8 1/4/ 
Iterating T', for large n 



T' 



cos7r/3, we find 
\ 



= 1/2 = 
/I 

= U 

V -1/8/ 
the off-diagonal element of 
the Jordan block starts dominating, hence (L|T'"|R) — > 
nUi,i(U-i)2,2 = 4n/9, so Dz{l/2) = 9/16. Similarly, 
we calculated for all m < 6 (note Dz{A) = Dz{—Aj): 



A_ 

7r/2 
7r/3 
7r/4 



Dz (cos A) 



A 



9/16 
1/3 



n/5 
2tt/5 
7r/6 



5(5- 
5(5- 



Dz {cos A) 



/5)/64 
-\/5)/64 
3/20 



(20) 



For A = 1, we have (L|T"|R) cx n^, so Dz\n oc n^^ -> 
0. For all other values of m < 300 we estimated ( 19 ) 



numerically and plotted the graph of Dz{/^) in Fig. [2| 
which clearly exhibits a fractal structure. This seems to 
be the first appearance of fractal behavior in an infinite 
strongly correlated system and calls for further analysis. 



Our numbers ( 20 1 have also been compared to the 
Bethe ansatz result [20, Dba = 1 - sin(2A)/(2A), for 
A ~ tt/iti, which overshoots the values of the table (20) 
for 4%, 9%, 13% and 15% for m = 3,4,5 and 6, respec- 
tively. However, since the assumptions used in Ref. [20] 
have been debated, it seems interesting to speculate 
whether our bound Dz may be saturating. 

Conclusion. We developed an original method to con- 
struct explicit matrix product ansatz for the weak cou- 
pling limit of NESS in the XXZ spin chain which is driven 
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FIG. 2. Mazur bound Dz on spin stiffness versus anisotropy 
A for an infinite XXZ chain, obtained by iterating (m 4- 1) x 
(m-|-l) transfer matrices at A = cos(7ri/m). Red dots indicate 
analytic results ( |20[ ) for m < 6. The inset magnifies Z)z(A) 
in a smaller range demonstrating its fractal nature. 



out of equilibrium by Lindblad baths attached at the ends 
of the chain. As a byproduct of this result we discovered 
a pseudolocal conserved quantity which - unlike previ- 
ously known local conservation laws - is not orthogonal 
to the spin current in the zero magnetization sector. Em- 
ploying this conserved quantity enabled us to rigorously 
estimate the spin Drude weight in the metallic regime 
|A| < 1 and prove the ballistic high-temperature trans- 
port there. Since the XXZ model can be mapped to a 
t-V model of interacting fermions, our result should be 
equally relevant for understanding the charge transport. 
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